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The Kugo-Ojima color confinement criterion, wliicli is based on the BRST sym- 
metry of the continuum QCD is numerically tested by the lattice Landau gauge 
simulation. We first discuss the Gribov copy problem and the BRST symmetry on 
the lattice. The lattice Landau gauge can be formulated with options of the gauge 
field definition, C/(link)-linear type or log (7 type. The Kugo-Ojima parameter ii^ 
which is expected to be — 1(5JJ in the continuum theory is found to be — 0.7<5J^ in the 
strong coupling region, and the magnitude is a little less in the weak coupling re- 
gion in log U type simulation. Those values are weakened even further in [/-linear 
type. The horizon function defined by Zwanziger is evaluated in both types of 
gauge field and compared. The horizon function in the log U version is larger than 
the other, but in the weak coupling region, the expectation value of the horizon 
function is suggested to be zero or negative. 

1 Introduction 

There are essentially two aspects in the manifestation of color confinement in 
the Landau gauge QCD. One aspect is the linear potential between quarks, 
which was conjectured by Gribov as a consequence of an enhancement of the 
singularity of the ghost propagatoci due to the restriction of the gauge field A 
on the transverse plane. Another aspect is the absence of free single colored 
particle state in the asymptotic Hilbert space, which culminates in the Kugo 
and Ojima color confinement criterion based on the BRST(Becchi-Rouet-Stora- 
Tyutin) symmetry: i.e. in the Landau gauge, a coefficient in the two-point 
function produced by the ghost, the antighost and the gauge field becomes 
—(5^, where a and h specify the color in the adjoint representation. Analytical 



calculation of this value is extremely difficult and so far no verification has been 
performed. In 1994, Zwanziger developed a lattice QCD theory for Gribov 
ambiguity. He claimed that, if the restriction to the fundamental modular 
region is achieved, the gluon propagator at momentum should vanish in the 
continuum limit. 

As the matter of first principle, it is still unclear whether the bases of 
these two theories are consistent or not, that is, 'Does the Kugo-Ojima the- 
ory properly resolve the problem of Gribov copy?' and/or 'Is the Zwanziger 
theory satisfactorily well defined in the continuum limit?' Considering this 
fundamental problem on one side, we make numerical verifications on conjec- 
tures of these theories on the other. In the lattice QCD test, we address the 
following problems: whether the gluon propagator is infrared finite, how sin- 
gular the ghost propagator_is in the infrared region, whether the Kugo-Ojima 
color confinement criterioila is satisfied, and how close to zero from negative 
the Zwanziger horizon function is. 

2 The Gribov problem and the lattice simulation of the Landau 
gauge QCD 

2. 1 The path integral formulation of the gauge fixed theory in the presence of 
the Gribov copy (Fujikawa, Hirschfeld) 

First we give a brief review of the palJi integral formulation of the gauge 
fixed theory by Fujikawa^ and HirschfeloB. They discussed a possible situation 
which may give a way out of the gauge fixing degeneracy problem (the Gribov 
problem). In the following, all indices representing multi-degree of freedom are 
suppressed. Let us define a gauge unfixed partition function as, 

Z = I dUe-^'^^'l (1) 

In derivation of the Faddeev-Popov(Faddeev-Popov) formula of the gauge 
f{U) = 0, one considers the determinant function A(C/^), 

^ ' V dg ) \ dA9 dg J 

where = Af^ {U), and one finds that A([/^) is indeed a function of 11^, and 
looks at the integral on the gauge orbit, , 

NiU) ^ [ dgA{U^)5{f{U^)). 
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Obviously N{U) is an orbit function, that is, N{U^) ~ N{U). At the intersec- 
tion points of the gauge orbit with the surface f{U) — 0, i.e., g — gi{U), 
Gribov's copies, the above delta function is transcribed to give 

N{U)= J d5^^|^%-5.(C/)) = J dgY,sign{A{mmg-g^,{U)). 

(2) 

These delta functions contribute or ±1, and thus, if the orbit function N{U) 
is non vanishing over all orbits, N{U) ^ 0, then in use of the identity 

1 = ^ / dgA{U^)d{fm), 

the standard FP procedure applies, and one factors out the gauge volume and 
obtains the formula for expectation values of functions, F{U), 

_ /dt7A(t/)^(/(t/))F([/)e-^^(^)/jV([/) 
{^{U)}\f„- jdUA{U)Sif{U))e-PSiu)/N(^u) ' ^' 

provided Fill) is gauge invariant. In the case when F{U) is gauge non- 
invariant, F{U) in (^) should be formally replaced with a gauge invariant 
function, 

j dg 

if one really wishes to obtain the Boltzmann average of the gauge non-invariant 
function F{U). The formula (|) allows the BRST formulation as in the fol- 
lowing, and if N{U) is particularly a constant, then the standard Lagrangian 
of the continuum theory is justified. The formula also derives a natural 
lattice simulation algorithm of the gauge fixed theory. The treatment of the 
gauge fixed theory here is, however, somewhat too formal for its practical use, 
particularly with respect to gauge non-invariant function in the presence of 
Gribov copies. We discuss this point in the following section. 

2.2 The simulation algorithm for the gauge fixed theory (Mandula-Ogilvie) 
and the Gribov problem 

Multiplying by the gauge volumes, the denominator and the numerator of the 
formula respectively, one recovers the Z in the denominator, and has in 
the numerator, 

JdU J dgA{U3)5{f{U3))F{U3)e-'^^^^^/N{U). (5) 
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As in (0), one obtains the numerator as 

dUY,si9n{A{U3^))F{U3^)e-P'^^^^/N{U) (6) 

i 

where U^^ is thciJLh Gribov copy on the orbit U^. This gives us the algorithm 
in the simulatioitJ 

where the last averaging () is that of simulation, i.e., that with respect to 
the Boltzmann weight e"'''^, and F{U) is a sign-weighted average of F{U^') 
on the gauge orbit U^. Readers should have already noticed that (0) gives 
trivial results when F(U) is gauge invariant, that is, we do not need any gauge 
transformation at all. Necessity of the function F{U) is rather impractical as 
well, and therefore we start from the formula (|^) with a gauge non-invariant 
function F(U) for work purpose. Then the result (0) is highly nontrivial, 
and F{U) can be viewed as being defined as a new gauge invariant function 
from Flu) and f{U). The sign- weighted average is impractical and a simple 
averag(o may be similarly considered as such. If the averaging has, however, 
a large fluctuation this formula would rather be interpreted as representing 
Gribov noise apart from the original function F{U), although could be 
considered as giving a sort of gauge invariant function out of F{U) and f{U). 
A modification of the gauge such that one chooses a unique copy among others 
on the orbit is favored, that is, a new gauge without Gribov copy, and 
in that case, the above formula is useful in practice, and this procedure may 
naturally change the result in the above formula (|^) without Gribov noise. 

2.3 The BRST formulation and the Gribov problem 

The obvious standard FP formula allows BRST formulation, 

{Fm\j „ - ^ 1 d^,eM5 J c/(t/))F({/)e-^^(^)/7V(C/), (8) 

Z' = J df^cMS J cf{U))e-^'^^^/N{U), (9) 
where 6 stands for BRST transformation, 

6c = iB, 6Af,{U) = Df,{U)c, 5c=-cc, 6B = 0, (10) 
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and the measure is defined as dfj, = dUdBdcdc, where dB integration is per- 
formed on R and dcdc is suitabUy defined differentiations with respect to Grass- 
mann numbers, c and c , such that 

dcdc 6^^^^= = det M, (11) 



dcdc c^Cye"""'" = (det M){~M-^)^y. (12) 

It was shown by Neuberger that if the gauge fixing function f{U) is a 
smooth function of compact variables ?7, then the expectation value of gauge 
invariant function F{U) becomes an indefinite fornid, 

{Fm\f gauge - ^, (13) 

which implies that all Gribov copies contribute to give total cancellation, in 
other words, the assumption that the above N{U) is non vanishing on all orbits, 
does not hold, and the formula (^ is totally meaningless. 

An essential point of his argument is as follows. Let us consider a general 
expression 

Z'{t) = / d^iex-p{tS I cf{U))G{U) (14) 



with a gauge invariant function G{U) as the Boltzmann weight. Then one 
finds that 

'^^'{t) ^ I d^x{5 [ cf{U))eMtS I cf{U))G{U) (15) 



dt ^ ^ 

can be written from nilpotency of 5, i.e., 6^ — 0, and from 5G{U) = as 
dZ' 



it) = J dfi6 (^( J cfiU))eMtS J cf{U))G{U) j . (16) 

If W{U, B,c,c) is an analytic function of the compact variables U, then one 
can show that 

dfidW{U,B,c,c) =0. (17) 

Thus it follows that Together with Z'(0) = / • ■ ■ dcG{U) = 0, one finds that 
Z'(t) = 0, so Z'{1) = 0. This concludes total cancellation of the Gribov 
copies of gauges given by analytic gauge functions of U. Thus one is forced to 
consider non-analytic gauge functions of U as desired gauge functions. As a one 
dimensional U{1) toy example which avoids (|l7|), one can consider U = e^^, 
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A = ImlogfJ, where A is not continuous at U = —1, and with definition 
6 A = c, one finds that 



//>7r /• oo 

dUdBdcdc 5{e-^^cA) = / dA I dBdcdc e^^'cc 7^ 0. 
J — 'K J —00 



(18) 



Although this problem is still open in the lattice BRST formulation, use of 
non compact variables in gauge fixing functions may be helpful. 

2.4 The Landau gauge and the Gribov problem 

Now we focus on the Landau gauge in SU{3) lattice QCD, that is. ,f(U) = 
dfiAfj_{U), where there are some options of definition Af^{U) as 

1. J7-linear OnJ; A^^,, = ^(C^x.p - U^,^)\traceless part, 



2. use of exponential majEl. U^^f^ = exp^^.^, A^. ^ — —A^^^^, where 
absolute values of all eigenvalues of A^^^j, do not exceed 47r/3. 

In the latter definition, Ax^p,{U) is not analytic with respect to compact variable 
Ux.fj. contrary to the^kirmei one. In both cases the Landau gauge, dA^ = 0, 
can be characterizecS'u in use of optimizing functions Fij{g) of such that 
5Fu{g) = for any bg. 

1. t/-linear definition; F^{g) - ^^y^^^ E.,^ (l " ^RetrC/^ J . 

2. use of exponential map; Fv{g) = \\A^\\^ = dV(N-^-i) E^.p tr {as.xJAI^^, 
where N is the number of colors and d is the dimension. It is noteworthy 
that 11^® IP is a continuous function of compact variables in spite of 
non analytic property of Ax.^i{U^^. 

In both options of the gauge field definition, the variation of the optimizing 
function, Fu{g), under infinitesimal gauge transformation g^^Sg — e, reads as 

t.Fu{g) = -2(9AS|e) + (e| - dD[A^)\^) + ■ • ■ , (19) 

where D^{A) denotes the covariant derivative in each definition; 

1. in ?7- linear version. 



D^iU,,,.)4> = - <^ ^i^i— 9^0 + [^,,^, 0], (20) 




where 

df,(f>^(f>{x + fi)-(l>{x), ^=^{cP{x + fi) + (f){x)), (21) 
2. in A = \ogU version, 

D^,{A,,,M = [^:r,M,<^], (22) 

where = adjA,,^ = and 

Gribosz. copy is generic phenomenon in both definitions as well as in the 
continuurm, there exist a lot of local minima of Fu{g) along the gauge orbit . 
Thus the naive Landau gauge loses its solid basis both in the theoretical and in 
the simulation view points for examination of gauge non-invariant quantities 
such as gluon propagator, ghost propagator, etc. 



2.5 Sophisticated Landau gauge 

Zwanziger devised various regions of the transverse plane of A^, i.e., dA^ = 0, 
depending on properties of a point on the plane. For example, one defines 
Gribov region D, as 

Vl = {A\- dD{A) > , = 0} , (24) 

where D^{A) denotes the covariant derivative, and —dD > implies that the 
Faddeev-Popov operator —dD is positive definite. A point on the Gribov region 
is a local minimum of Fu{g), but it is known that some points on the Gribov 
region can be gauge copies of each other. Thus one defines the fundamental 
modular region A as the absolute minimum along the gauge orbits. 

A= {A\\\A\\^ ^Ming\\As\\^}, A c n . (25) 

Then one defines a corresponding region Ua of configuration U , as 

Ua = {U\A{uf e A} 

{A^ denotes transverse component of A(U)) (26) 
Putting an indicator function 9a of the set Ua as 

6a{U) ^lif U eUA, 9a{U) = Oif U (^Ua, 
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and putting Oa{U) 
as 



= 1 — 9a{U), one can define a corresponding gauge function 
fAiU) = dA-9A{U)+h{U) (27) 



Tlie arguments in tfie preceding subsections formally applies for a correspond- 
ing gauge function containing non-analytic Heaviside function. The gauge- 
fixing algorithm in the simulation is required to attain the absolute minimum 
of the Fij(g) along the gauge orbit. But the global minimization is difficult in 
general and developing the efficient algorithm is still an open problenxHI. 

3 The Kugo-Ojima confinement criterion and the Gribov-Zwanziger's 
theory 

3.1 Kugo-Ojima' s theory 

A sufficient condition of the color confinement given by Kugo and OjimJi is 
that u"'' defined by the two-point function of the FP ghost fields, c(a;),c(j/), 
and A^{y), 

J e'P^-'y\0\TD^c''{x)g{A, x cfiyMd^x = {g^, - Pl^)u''\p^) (28) 
satisfies w"''(0) = 

Brief survey of Kugo-Ojima's argument that = —5^ is a sufficient con- 
dition of the colour confinement is the following. 

1. The BRS transformation is given by replacing gauge transformation 
parameters 0° as 9°'{x) —>■ Xc°'{x), where A is an imaginary Grassmann 
number and c"'{x) is the ghost field. Putting C{x) = gc"'{x)Aa where 
Aq is antihermitian such that [Aa,Af,] = fabc^c and Sb = XSb, and 
antihermitian matrix Afj_(x) — gA1^{x)Aa, the BRS transformation reads 

SB(b - SbA^ - D^{A)C - d^C + [A^,C]. 

The nilpotency requirement of Sb derives, 

For each c°'{x), one introduces an anti-ghost c°(a;), and similarly the 
matrix C, and defines, 

SbC = iB, 
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then the nilpotency gives 6b B = 0. Then the gauge- fixing and Faddeev- 
Popov Lagrangian can be written as 

jCg:F+:FV = -^<5B[c«(9M« + ^aB^] 

= B^d^^Al + iaS'^B'^ + i^d^D^c". (29) 

Together with the original gauge invariant Lagrangian, the BRS invari- 
ance of the total Lagrangian follows from the nilpotency 5b^ = 0. Thus 
the BRS symmetry is a symmetry of a gauge fixed Lagrangian, so to say, 
a quantum gauge symmetry. The corresponding conserved BRS charge 
reads as 

B^Doc" - ^oB" ■ + '^-gdo^ ■ (c x c)" 

where (F x G)" = fabcF^G'^- Thus one defines physical space VpUys = 
{\phys)} as 

Qslphys) = 0. 

It is to be noted that both ghost and anti-ghost fields are considered as 
hermitian fermi fields so that the derived Hamiltonian should be hermi- 
tian. This assignment is important for derivation of the unitarity, but 
due to this assignment, there necessarily involves an indefinite metric in 
a space of ghost and anti-ghost fields. 

2. The BRS algebra is given by BRS charge, Qb, and FP ghost charge, Qc, 
as 

Ql = 0, [iQc, Qb] = Qb, [Qc, Qc] = 0. 

Since these Qb,Qc are commuting with other conserved charge, all asymp- 
totic one particle states can be classified by irreducible representations. 
Due to the nilpotency of BRS charge, Qb, there only exist BRS singles 
and doublets. From hermiticity of ghosts, Qc is defined as a generator of 
scale transformation of FP ghosts. It is to be noted that while Qc being 
hermite, FP ghost number is counted by Npp = iQc, and that among 
FP ghost number eigenstates, {M\N) ^ only if M = —N. Due to this 
metric structure, BRS doublets always appear in pair of opposite sign FP 
ghost numbers, and this pair is called a BRS 4-tet (quartet). Under the 
assumption that BRS singlets have positive metric, it is proved that 
Vphys has positive semidefinite in such a way that BRS 4-tet particles 
appear only in zero norm. 



Ql 



(Px 
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3. Prom the Ward-Takahashi identities, 



F.T.mTD^c'^{x)c\ym = -*(0|rA"(x)i?^(y)|0)] 



(30) 



where F.T. imphes the Fourier transformation, it follows that Heisenbcrg 
operators, iP^c", c", i?" necessarily have massless asymptotic fields 
when xq — > ±00, 

B'^{x)^(3\x) + ---, 
D^c'^ix) ^ 9^7"(a;) + ■ • • , 

if{x)^r{x) + ---- 

One finds from the BRS transformation that for each colour a, a set of 
the above massless asymptotic fields form a BRS 4-tet (quartet) such 
that 

[iQb,x''{^)]=i''{x), {iQB,nx)}=if3{x) 
{iQb,i''{x)} = 0, [iQb,I3''{x)]=0 

and 

mx),^'{y)} = -i[r{x)J'{y)] = -S'^'Dix - y) 

4. With respect to symmetry breaking in general, the following statements 
are equivalent; 

(a) Charge of symmetry, Q = J d^xjo, is well-defined. 

(b) Symmetry is not broken, Q\0) =0 

(c) There exist no massless one particle states in spectrum, 

{0\jf.{x)\p' =0)=0 

5. The Noethcr current corresponding to the conservation of the colour 
symmetry is gJ^ — d'^ F^^ + {Qb^D^c}, whore its ambiguity by diver- 
gence of antisymmetric tensor should be understood, and this ambiguity 
is utilised so that massless contribution may be eliminated for the charge, 
Q", to be well defined. 
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6. Denoting g{A^ x c)" u'^d^'j'', and then D^c" ^ (1 + u)'^d^j'', one 
obtains that 

F.T.{0\TD,c'^{x)9{A^ x c)''(y)|0) = {g,. - ^Kip'), 

provided has a vanishing expectation value. The current {Qb,D^c\ 
contains the massless component, (1 + u)1^d^(3^{x). We can modify the 
Noether current for colour charge such that 

In the case of 1 + u = 0, massless component in gj'^ is vanishing and the 
colour charge 

Q- = j d?x{QB.g-^Do(^{x)} (31) 
becomes well defined. 

7. The physical state condition Qs^phys = together with the equation 
(|l|) implies that all BRS singlet one particle states |/) S Vphys are 
colour singlet states. This statement implies that all coloured particles 
in Vphys belong to BRS 4-tet and have zero norm. This is the colour 
confinement. 

8. In the course of their derivation, they assume Lorentz invariance and 
that the colour symmetry is not broken. 

9. They also proved that if the vector massless asymptotic field is missing 
in a channel a, and if the channel a belongs to the image of 1 + m then the 
massless 4-tet in j° can not be cancelled, and the colour symmetry with 
charge Q", is spontaneously broken. (Inverse Higgs mechanism theorem) 

The corresponding Euclidian expression is as follows. 



which can be calculated by 




ir^\^^D,^[-AA'Vj j^iSp.-^K\p% (32) 

where A" is a normalized antihermitian basis of Lie algebra, V a lattice volume, 
and the ghost propagator is given by 



(^^^~<^ = {''{^''^>^'] }■ (33) 

xyl 



dD 
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3.2 Zwanziger's theory 

The fundamental modular region A is specified by the absolute minimum 
along the gauge orbits in the Gribov region fl. 

A = {A\\\Af =Mmg\\Aof}, A c Q = {A\ - dD > , OA = 0} . (34) 

Zwanziger relaxes the periodicity restriction on the gauge transformation 
g, and imposes larger periods than the original. Then some two points in 
the fundamental modular region A may be bridged to be Gribov copies of 
each other, and one of them is not the absolute minimum of the minimizing 
function along the gauge orbit anymore. Surviving points as the absolute 
minimum consist of core region 5 ( S C A). In the so defined core region S, 
a horizon function H{U) given below is negative. 

The Horizon function is defined as follows. Let two point tensor G^^xy"''' 

be 

G^uxv''' = tr[x^^D^^{-D,)X^^ . (35) 
Then H{U) is given as 

HiU) = J2 G.^xy'"' - {N' - l)EiU) (36) 

x,y,a 

where E{U) reads as follows; 

1. in J7-linear version, E{U) = J^i ^^Ui, 

2. mA = logC/ version, E{U) = ^ ^tr (^X"-^ S{Ai)X''y 



where Ai = adjAi , and S{x) 



l,a 

x/2 



th{x/2) ' 



Let us define an average tensor Gf^,^xy be G^i^xy^""^ = provided color 

symmetry is not broken. One sees that a Fourier transform of the average 
tensor, 

(7 — p-ip{x-y)Q 



V 

x,y 



^ IJyVxy 



takes a form 

G^.W' = (I) ^^"^ - [S,. - ^) (37) 
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where e = {E{U)) /V , and that it is related with the horizon function as 



lim Gf,f,{p) 



(38) 



He defined the augmented core region 'f = {U : H{U) < 0} n il 
(Sc^'Cfi). 'J and A are qualitatively similar, and he defined the 
partition function in the path integraLin use of the corresponding Lan- 
dau gauge function and concludecEl in the infinite volume limit that 
limy^oo {H{U))/V — 0. Putting Kugo-Ojima parameter as u°-^{0) = — (5°^c, 
one finds from (37), (|3|), that 



Q+(.-l)c-e=(.-l)(c-^)=0, 

which is called horizon condition. Since we can measure c and e by the lattice 
simulation, we can check to what extent Zwanziger's horizon condition holds in 
our simulation. With respect to the value e/d, note that the classical vacuum 
is characterized by e/d = 1. 



4 The numerical results 

4-1 The algorithms of lattice Landau gauge fixing 

In the definition A — losi/, our Landau gauge fixing algorithm is as follows. 
We define the gauge fielcQ on links as an element of SU (3) Lie algebra as, 

e^-'" = U^,^ , where A^^ = -A^.,^. (39) 

We perform the gauge transformation as e"*^ *' = gt^^''''^ 9x+fi and define 
\dA\ = Max^.^JdAll \\dA\\^ = ^ i^a^taA,. The Landau gauge 

is realized by minimizing 1 1 | p via a gauge transformation g''Ug, where g = e"^ . 
In order to obtain e, we switch the following two methods, depending on the 
current value of \dA\ in comparison to some critical parameter \dA\cr- 

1. When \dA\ > \dA\cr, £x = ng^n dA^ with suitable parameter 1 < 77' < 2.2 

2. When \dA\ < \dA\cr, e = i-d^D^{A))-h]dA where 1 < < 2 is a 
parameter. 

The restriction to the fundamental modular region is not always achieved. 
But, we observed that the obtained ncpp \\A\\ is larger or smaller than that 
obtained after the smeared gauge fixingJ within 1% accuracy. 
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In case of tim. {/-linear definition of gauge field A, we perform the site-local 
exact algorithnxJ, with suitable over-relaxation parameter, rj = 1.6, starting 
from gauge fixed configurations oi A = log U. It is found that the exact over- 
relaxation g = W^'^ is faster than the stochastic over-relaxation, where W 
is the site-local exact solution obtained by solving the nonlinear equation for 
finding the best fit gauge transformation on even-(odd-)site. 



^.2 The Kugo-Ojima two-point function and the ghost propagator 
The FP operator is 

M[U] = -id-D{A)) = -iD{A)-d), (40) 

and we define Ad{Af,) by putting D^{A) = di_,+Ad{Af,). The inverse, M^'^[U] = 
{Mq — Mi[U])^^, is calculated perturbatively by using the Green function of 
the Poisson equation Mq^ = {~d'^)~^ and Mi = df_iAd{A^) , as 

= Mo' + J2 {Mo'MifMo^- (41) 

k=0 

The ghost propagator ( p3[ ) is infrared divergent and its singularity can be 
parameterized as p"^'^^+"\ where p^ = ^(4sin^ < < L/2). 

It depends on (3 slightly, but its finite-size effect is smal0. These qualitatiiie 
features are in agreement with the analysis of the Dyson-Schwinger equatioit-3. 




1 2 3 4 11.522.533.54 



Figure 1: The ghost propagator as func- 
tion of the lattice momentum. The data are 
/3 = 5.5(box) and 6.0(triangle), 16*. The 
fitted curve is 1.287/p2-779 f^j, ^ ^ 5 5 
1.162/p2-545 ^ ^ g g (dashed). 



Figure 2: The Kugo-Ojima parameter |tij^| 
as the function of the spatial extent of the 
lattice aL(/m). The data are /3 = 6 (tri- 
angle), and /3 = 5.5 (box) 8"*, 12-4, le* from 
left to right, respectively. 
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We measured the lattice version of |m^(0)I on 8'^ x 16,8* 12"* and 16* for 
P = 5.5 and 6. When f3 — 6 and the lattice size is small, the Polyakov 
loop distribution deviates from the uniform distribution. In this case, we 
perform the Z3 rotation by multiplying the global phase e^"r' such that the 
distribution concentrate around one angle, before we measure the Kugo-Ojima 
two-point function. 

We obtained that ug(0) is consistent with —cS^, c — 0.7 in SU{3) quenched 
simulation, /3 = 5.5, on 8*, 12* and 16*. 

4-. 3 The gluon propagator 

The gluon propagator is infrared finite. We parameterized the zero-temperature 
lattice data using the Stingl's Factorised Denominator Rational Approximant 
(FDRA) method. The effective mass of the gluon in the analysis of 8^ x 16 is 
found to be about 60QMeV. 

The infrared finiteness is in accordance with the Kugo-Ojima color con- 
finement mechanism. As stated in the their inverse Higgs mechanism theorem, 
if we have no massless vector poles in all channels of the gauge field- and 
if the color symmetry is not broken at all, it follows that 1 -|- m = Ota. 

5 Discussion and conclusions 

We performed the first test of the Kugo-Ojima color confinement criterion in 
lattice Landau gauge. We observed that the 8*, 12* and 16* lattice data. The 
data of /3 = 5.5 indicates that u = —0.7, and those of /? = 6.0 are smaller by 
about 10%. 

In the Zwanziger's theory, the two-point function Gui/(fc^) can be expressed 
in terms of the Kugo-Ojima two-point function as (|37|). Zwanziger's horizon 
conditioril in the infinite volume limit reads as G'p^(O) = e = (-^^7^). In terms 
of the Kugo-Ojima parameter c, the left hand side can be written as (e/4) + 3c, 
and the horizon condition is that c = e/d. In the table below ei and 62 stand 
for e in our 16* lattice simulation of the first and the second option of the 
gauge fields, respectively. 

If the gauge fixing could be performed so that it brings the configuration 
into the core region or the augmented core region and if the infinite volume 
limit is considered somehow, then the legitimate check of the horizon condition 
could be done. The core gauge fixing is, however, difficult, and even impossi- 
ble in general, which implies that the core gauge is literarily not the gauge. A 
configuration of the core region of period L, belongs to a fundamental mod- 
ular region of larger period NL as well as to a fundamental modular region 
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of period L by definitiom. Such a configuration is particular one in the fun- 
damental modular region of period L, to which a generic configuration in the 
fundamental modular region of period L cannot be gauge transformed even 
by a relaxed gauge transformation of larger period. A relaxed gauge transfor- 
mation with larger period NL unique up to global gauge transformation can 
bring the generic configuration above to the fundamental modular region of 
period NL, but at the same time breaks the periodicity of L in general. Thus 
restriction to core region is neihter an argument of 'gauge' nor an argument 
of trace disappearance of periodicity. It is highly dynamical hypothesis that 
the core region and the fundamental modular region give the same limiting 
correlation fucntions, i.e., dynamics, in the infinite volume limit. Thus we take 
a standpoint that the horizon condition derived from Zwanziger's restriction to 
the core region is simply tested by the dynamics of the fundamental modular 
region, and give the direct results in the table although obviously not in the 
infinite volume limit. 

Table 1: /3 dependence of the Kugo-Ojima parameter c, the tensor G^fi/A, trace e divided 
by the dimension d. The suffix 1 corresponds to the n-Unear and 2 corresponds to the log U 
definition. Data are those of 16^, except /3 = 5.5 (7-linear data, which are those of 8*. 





Cl 


ei/d 


G^^i(0)/d 


C2 


62 /d 


G^M2(0)/d 


5.5 
6.0 


0.570(58) 
0.576(79) 


0.780(3) 
0.860(1) 


0.622(45) 
0.647(57) 


0.712(18) 
0.628(94) 


0.657(1) 
0.693(1) 


0.698(14) 
0.644(70) 



Simulation data show in general that when (3 becomes larger, e becomes 
larger, while c has an opposite tendency. This fact itself does not necessarily 
disprove the horizon condition, but our data of c which is calculated in the 
A — \ogU version already gives the zero-intersection of G^^/d — 62 /d in the 
increase of [3 from 5.5 to 6. 

This work was supported by KEK Supercomputer Project (Project No.OO- 
57), and by JSPS, Grant-in-aid for Scientific Research(C) (No.ll640251). 
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